I 



LOAN COPY: RETURN TO 
AFWL (WLIL-2) 
KIRTLAND AFB, N MEX 


BACK FLOW FROM 
JET PLUMES IN VACUUM 


by Norman T. Grier 

Lewis Research Center 
Cleveland, Ohio 





NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


WASHINGTON, D. C. • JANUARY 1969 


0X31550 


TECH LIBRARY KAFB, NM 

nniniiin 

□131558 


BACK FLOW FROM JET PLUMES IN VACUUM 

By Norman T. Grier 

Lewis Research Center 
Cleveland, Ohio 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 2215 1 — CFSTI price $3.00 


V 



ABSTRACT 


Approximate calculations were performed to determine the increase in particle 
density at points surrounding a jet nozzle due to the exhausted gas in the jet plume. 

The calculations were made for helium, argon, and nitrogen for a nozzle throat radius 
of 1 centimeter. The chamber density ranged from 10“® to 1 gram per cubic centi- 
meter. The plume flow is assumed radial with the source point at the nozzle exit on 
the nozzle axis. 


BACK FLOW FROM JET PLUMES IN VACUUM 
by Norman T. Grier 
Lewis Research Center 

SUMMARY 

Approximate calculations were performed to determine the increase in particle 

density at points surrounding a jet nozzle due to the exhausted gas in the jet plume. 

The calculations were made for nitrogen, helium, and argon for a nozzle throat radius 
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of 1 centimeter. The chamber density ranged from 10 to 1 gram per cubic centi- 
meter. The plume flow is assumed radial with the source point at the nozzle exit on the 

nozzle axis. The results indicate that for a chamber density of 1 gram per cubic centi- 
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meter, the back flow density could be as high as 10 to 10 particles per cubic centi- 
meter. 


INTRODUCTION 

For spacecraft accelerating or maneuvering in space, the back flow from the jet 
exhaust may be a large contributor to the ambient environment of the spaceship. If this 
exhaust gas is detrimental in any manner to the planned purpose of the flight, special 
shielding may be required. In order to get an estimate of the density of the back flow, 
an approximate calculation is performed in this report. The gases used are helium, 
nitrogen, and argon. In all the calculations the jet exit Mach number M g is 5. 

The results presented in this report may be considered an extension of a companion 
paper published in the Proceedings of the Sixth International Symposium on Rarefied 
Gasdynamics (ref. 1). In that paper, results for points located up to 100 exit radii in the 
plane of the nozzle exit are presented. In this report, results are presented for points 
located up to 1000 exit radii in the plane of the nozzle exit as well as for points in planes 
ahead of the nozzle exit. Also, in the previous paper (ref. 1), the maximum angle that 
the plume flow could turn through was assumed to be the Prandtl-Meyer expansion angle. 
This assumption is not made herein. Rather it is assumed that the analysis of Hill 
and Draper (ref. 2) is valid throughout the plume flow field. And finally, the analysis 
presented herein is more detailed than that of the previous paper (ref. 1). 



The gas in the plume is considered to be continuum until the density is such that the 
mean free path is comparable to the distance the particular mass of gas has travelled 
from the nozzle exit. That is, the jet flow is continuum until the Knudsen number has a 
value near one. Afterwards, the flow is assumed free molecular with no interparticle 
collisions occurring. This model is shown in figure 1. 

Since most of the constant Knudsen surface is usually far from the jet exit plane, it 
is assumed that the flow has nearly reached its adiabatic limit and has a radial direction, 
that is, the streamlines diverge from a common source point located at the nozzle exit. 
The validity of these assumptions is discussed in references 2, 3, and 4. The calcula- 
tions were performed for a nozzle throat radius of 1 centimeter, a value of Knudsen 
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number of one, and chamber densities ranging from 10 to 1 gram per centimeter. For 
a chamber temperature of 300 K, these densities correspond to chamber pressures for 
helium of —6x10”^ to 6x10^ atmospheres, for argon of — 6xl0~® to 6x10^ atmospheres, 
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and for nitrogen of —9x10 to —9x10 atmospheres. In all the calculations the environ- 
ment excluding the exhaust gas was assumed to be a total vacuum. A high-speed digital 
computer was used for the calculation. 


ANALYSIS 

As mentioned earlier, references 2, 3, and 4 show that far from the nozzle exit, 
the exhaust from jet nozzles in vacuum can be closely approximated by radial flow. For 
radial flow, the mass flux pV varies inversely as the square of the distance from the 
source point. Far from the source, the flow has nearly attained its adiabatic velocity 
limit, so the density decrease along each streamline is proportional to the inverse 
square of the distance from the source. Applying this description to jet plumes, refer- 
ences 2, 3, and 4 show that the density along the plume axis ( 9 = 0) is approximately 
given by 



where r* is the nozzle throat radius, r is the distance from the source point which is 
assumed located on the nozzle axis in the plane of the nozzle exit, p c is the chamber 
density, and B is a constant that depends on the type of gas. (See fig. 2 for coordinate 
system. ) (All symbols are defined in the appendix. ) For r constant the density may be 
written approximately as a function of 6 as (ref. 2) 
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( 2 ) 


p = p e=0 e ' x2(1 - cos e) 2 


= 4p b/^V— \ e^ 2 * 1 - 008 e > 2 
C \ e/\ r / 


where r is the nozzle exit radius 
e 



( 3 ) 


B=-^-/x^y / Y-i-y /<y ■ l, ( 4) 

4C^V* V V'+V 

where C ^ is the thrust coefficient and C is the ratio of the velocity to the maximum 
velocity. The ratio, C, is assumed to be approximately 1 on the constant Knudsen sur- 
face. The thrust coefficient Cp was evaluated using equation (4. 33) of reference 5. 

Equations (1) to (4) apply as long as the flow is continuum. Actually a transition 
from continuum to free molecular flow occurs over a region in the plume. It is assumed, 
however, that this region can be approximated by a surface on which the Knudsen num- 
ber is constant. For the present study, the Knudsen number is defined by 

K =— (5) 


where l is the mean free path and is given in reference 6 approximately as 


l 


3. 736X10 


-25 M 



cm 


( 6 ) 


where M is the molecular weight in gram per gram -mole and 6 is the molecular diam- 
eter of the molecules in centimeters. Substituting equation (6) into equation (5) and 
solving for p yield 
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p = 3. 736x10 


M 


(7) 



Equating equations (7) and (2) results in the following equation of the constant Knudsen 
surface 


r 25 ® 

J_ = 1 . 07x10^° — - 

M 


— '\ 1//2 e - x2 ( 1 " cos e ) 2 

. A e/ 


( 8 ) 


In equation (8) the maximum value for 0 was taken as the value corresponding to a 
value of r/r of 0.1. Solving equation (8) for „„ gives 

0 IXlcLX 


r 




max 



(9) 


To determine the effect of limiting r/r g to values greater than 0. 1, a calculation 
was made using a L„ v corresponding to r/r. of 0. 05. The difference could not be 
discerned on the figures. Where there was any difference, it occurred for y/r g less 
than 20. 

On the constant Knudsen surface, the gas is assumed Maxwellian with a mean 
velocity V in the radial direction. From Noller's report (ref. 7), the density at a 
point P(x,y,z) due to the molecules leaving the elemental volume dr adjacent to the 
constant Knudsen surface is 


dn-r, = n ir 


-3/2, 


r u2 Jlu 


COS if/ + 


U 2 cos 2 i//je U cos + erf(U cos \p) 


( 10 ) 


where n T is the number density of molecules in dr, U - | v|/(2RT/M)^^ 2 , that is 
| V J nondimensionalized with respect to the most probable thermal velocity, i// is the 
angle between the direction of V and the direction of flight of the molecules contributing 
to the density at P(x, y, z), erf (U cos if/) is the error function, and d!2 is the element 
of solid angle of dr as seen from point P(x, y, z). 

Using the coordinate system shown in figure 1 gives 
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cos \p = L* 

|L| 


(ID 


dfi = ? • 11 dcr (12) 

|E | 3 

where L is the vector from dT to the point P(x, y, z), r is the unit vector in the radial 
direction, n is the unit vector normal to the constant Knudsen surface at dT, and da 
is an element of area on the constant Knudsen surface at dr. The coordinates for dT 
are denoted by (r, 9, <p) and the coordinates for point P by (x, y, z). Using figure 2 
and restricting x to be zero (the plume is symmetrical about the z-axis) result in 


— = r y X- sin 9 cos (p - — + — cos 9 j+ 9 ( — cos 9 cos cp - — sin 9 J - cp -3L sin (p 


y. 

r 


(13) 

where r, Q, and <p are the spherical coordinate unit vectors. Substituting equation (13) 
in equation (11) yields 


2— sin 9 cos cp - — + — cos 9 

COS rp = ;; (14) 

I L | 


Now on the constant Knudsen surface, r is related to 9 through equation (8), which can 
be written as 


where 


r = 


r fl=O e 


-X 2 (l-cos 9) 2 


T 9=0 


1.07X10 25 


Bp c r*K6 2 

M 



(15) 


(16) 
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Equations (15) and (16) define the constant Knudsen surface in terms of 0 as an 
independent variable. The other independent variable is the azimuthal angle <p. There- 
fore, (ref. 8) 


n do = — x — d0 d<p 
30 3 cp 


Now 



(1-cos 0)^ 


(17) 


(18) 


So that 


3f 

30 


r 0=O e 


-X 2 (l- 


cos 0)' 


[- 


r2X 2 (l 


cos 0) sin 0+0 


(19) 


= r = tpr sin 0 (20) 

3 (p 3 <p 


Therefore, 

n da = |^02X 2 (1 - cos 0) sin 0 + r J r 2 sin 0 d0 d <p 
Substituting equations (13) and (21) into equation (12) gives 


( 21 ) 


drt = 


i l 


sin 6 cos <p[ 1 + 2X* cos 6 - 2X^ cos^ 6 
r 


r 

r 


f( c 


+ _£_(cos 0 - 2A 2 sin 2 0 + 2A 2 sin 2 0 cos 


sin e d0 d( P 


( 22 ) 


The density at point P(0, y, z) due to molecules leaving the total surface is 
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n p = 7r"^ 2 e' 


2 r ’ /2 

u I / df2n ri 

• 4=0 J cp=-ir/2 


- U cos \}/ 
2 


+(i + oW*)e^Vl£j[i 


1 + erf(U cos ^/) 


(23) 


where L„, is given by equation (9), is given by equation (22), and cos ^ is given 
by equation (14). 

In order to get the results presented, it was assumed that the velocity V was 
95 percent of its adiabatic limiting velocity. This precludes the necessity of assuming 
a value for the temperature. The gas is assumed to be ideal. Therefore, the equation 
for the dimensionless velocity U is 



= 0.95 


y 


0.0975(y - 1) 

The molecular diameter 6 required for the various gases are (ref. 6) as follows: 


The number density n T is 


Nitrogen: 

3. 76xl0" 8 cm 

Helium: 

2. 18xl0 -8 cm 

Argon: 

3.67X10 -8 cm 


N o 

n=—p 
7 M T 


where N Q is Avogadro’s number and p T is given by equation (2) or (7) evaluated at 
dr. 
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RESULTS AND DISCUSSION 


A typical shape of the constant Knudsen surface is shown in figure 1. This figure 
is drawn to scale relative to the exit radius shown for helium using K = 1, r* = 1 centi- 
meter, and p c = 10~ 7 gram per cubic centimeter. This Knudsen surface extends to 
approximately 34 r e . If the chamber density had been 1 gram per cubic centimeter, the 
surface would be 10^ times larger. This follows from equation (16) where it can be 
seen that r^_Q varies linearly with the chamber density. 

In figures 3 to 5 is shown the back flow number density for helium, argon, and 

nitrogen, respectively, for values of y/r up to 1000. The results are given for 
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chamber densities from 10” to 1 gram per cubic centimeter, a Knudsen number of 1, 
and a nozzle throat radius of 1 centimeter. Because of approximations made in the 
analysis, the results are only expected to indicate trends and to give an order of magni- 
tude of the back flow density. The results for helium and argon show the effect of 
changing the molecular weight by a factor of 10. From figures 3 and 4 it is seen that in 
no case does the back flow density decrease by a factor of ten. In most cases the de- 
crease in changing from helium to argon is only by a factor of approximately four. 

The calculation for nitrogen shows the effect of changing from a monatomic to a 
diatomic gas. Comparing figures 3 and 4 with figure 5 shows that the back flow for 
nitrogen decreases much faster with y/r g than does helium or argon. Also it is noticed 
that the back flow density for nitrogen also decreases faster with chamber density than it 
does for helium or argon. As was indicated previously, the difference in molecular 
weights cannot account for this variation. So most of the difference in the behavior of 
nitrogen from that of helium and argon must be due to the difference in the value of y. 
This causes X and U to be different. 

One final point to be noticed from the figures is that the back flow in the plane of the 
nozzle exit continues to decrease as the distance from the nozzle y/r g increases. The 
back flow for points located in planes ahead of the nozzle exit, however, first increases, 
peaks, then decreases as y/r e increases. This can be seen from figures 3 to 5. The 
peaks are broadest for planes located the farthest from the nozzle exit. The calculations 
were made for planes located 0, 2, 5, and 10 exit radii ahead of the nozzle exit. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 23, 1968, 

124-09-19-01-22. 
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APPENDIX - SYMBOLS 


A* 

2 

nozzle throat area, cm 

n p 

particle density at point 

A e 

2 

nozzle exit area, cm 


P(x, y , z) , particles/ cm 2 

B 

defined by eq. (4), dimension- 

n T 

particle density in volume 

C 

less 

ratio of velocity to maximum 

R 

element dr 
gas constant 


velocity 


nozzle throat radius, cm 

C F 

nozzle thrust coefficient 

r e 

nozzle exit radius 


maximum nozzle thrust coef- 

r, 6,cp 

spherical coordinates 

max 

ficient 

r, e,(p 

unit vectors of r, 9, cp 

K 

Knudsen number, 1 / r 

u 

dimensionless velocity, 

L 

vector distance from dT on 


| V|/ (2RT/M ) 1//2 


constant Knudsen surface to 

V 

mean velocity of particles in 

l 

point P(x,y,z), cm 
mean free path of gas par- 

7 

jet, cm/sec 
isentropic exponent 

M 

ticles, cm 

molecular weight of gas 

6 

molecular diameter of molecule 
cm 

M e 

nozzle exit Mach number 

A 

defined by eq. (3) 

N o 

Avogadro's number, par- 

P 

3 

mass density, g/cm 

n 

ticles/g-mole 

particle density, particles/ 
3 

P C 
* 

3 

chamber mass density, g/cm 
angle between direction of V 

n 

cm 

unit vector normal to constant 
Knudsen surface 

and direction of L 


9 



REFERENCES 


1. Grier, Norman T. : Back Flow from Jet Plumes in Vacuum. Presented at the 6th 

Rarefied Gas Dynamics Symposium, Cambridge, Mass. , July 22-26, 1968. 

2. Hill, Jacques A. F. ; and Draper, James Stark: Analytical Approximation for the 

Flow from a Nozzle into a Vacuum. J. Spacecraft and Rockets, vol. 3, no. 10, 

Oct. 1966, pp. 1552-1554. 

3. Sherman, F. S. : Self -Similar Development of Inviscid Hypersonic Free -Jet Flows. 

Rep. 6-90-63-61, Lockheed Missiles and Space Co. , May 23, 1963. (Available 
from DDC as AD-618954). 

4. Albini, Frank A. : Approximate Computation of Underexpanded Jet Structure. AIAA 

J. , vol. 3, no. 8, Aug. 1965, pp. 1535-1537. 

5. Shapiro, Ascher H. : The Dynamics and Thermodynamics of Compressible Fluid Flow. 

Vol. 1. The Ronald Press Co. , 1953. 

6. Santeler, Donald J. ; Holkeboer, David H. ; Jones, Donald W. ; and Pagano, Frank: 

Vacuum Technology and Space Simulation. NASA SP-105, 1966. 

7. Noller, H. G. : Approximate Calculation of Expansion of Gas from Nozzles into High 

Vacuum. J. Vacuum Sci. Tech. , vol. 3, no. 4, July/Aug. 1966, pp. 202-207. 

8. Phillips, H. B. : Vector Analysis. John Wiley & Sons, Inc., 1933. 


10 






Back flow particle density, rip, particles/cm 



(c) For points located in plane 5r e ahead of nozzle exit (d) For points located in plane 10r e ahead of nozzle exit 

(i. e., z/r e = -5). (i. e. f z/r e = -10). 

Figure 3. - Concluded. 
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Back flow particle density, rin, particles/cm 





